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ABSTRACT 

We consider a model of axisymmetric neutron star magnetosphere. In our approach, 
the current density in the region of open field lines is constant and the returning cur- 
rent flows in a narrow layer along the separatrix. In this case, the stream equation 
describing the magnetic field structure is linear both in the open and closed regions, 
the main problem lying in matching the solutions along the separatrix (Okamoto 1974; 
Lyubarskii 1990). We demonstrate that it is the stability condition on the separatrix 
that allows us to obtain a unique solution of the problem. In particular, the zero point 
of magnetic field is shown to locate near the light cylinder. Moreover, the hypothe- 
sis of the existence of the nonlinear Ohm's Law (Beskin, Gurevich & Istomin 1983) 
connecting the potential drop in the pair creation region and the longitudinal electric 
current flowing in the magnetosphere is confirmed. 
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1 INTRODUCTION 

A few years after the discovery of radio pulsars (Hewish et al 
1968) the basic properties of the neutron star magnetosphere 
were actually clarified. Firstly, the importance of one-photon 
electron-positron pair creation in a superstrong magnetic 
field was shown (Sturrock 1971; Ruderman & Sutherland 
1975). It means that the magnetosphere of a neutron star 
must be filled with a plasma screening the longitudinal elec- 
tric field (Goldreich & Julian 1969), the charge density pe 
being nonzero. This screening results in the corotation of a 
plasma with the star. Clearly, such a corotation is impossible 
outside the light cylinder Rl = c/Q, where is the angu- 
lar velocity of a star. Hence, in the magnetosphere there 
form two essentially different groups of field lines: closed, 
i.e. those returning to the star surface, and open, i.e. those 
crossing the light cylinder and going to infinity. As a result, 
along open field lines the plasma may leave the neutron star 
and escape from the magnetosphere. The charge pe going to- 
gether with plasma produces a longitudinal electric current. 
It was also shown that it is the ponderomotive action of the 
electric current fiowing along open field lines and closing on 
the star surface that diminishes the angular velocity of a 
pulsar. 

At the same time, the equation describing the struc- 
ture of an axisymmetric stationary magnetosphere in the 
force-free approximation was obtained by several authors 
(Mestel 1973; Scharlemann & Wagoner 1973; Michel 1973ab; 
Okamoto 1974) . This second-order elliptical equation on the 



magnetic fiux function has the Grad-Shafranov form i.e. con- 
tains two 'integrals of motion', a longitudinal electric cur- 
rent and an angular velocity of magnetic surfaces, which are 
constant along magnetic field lines. Moreover, a few exact 
solutions of the stream equation were obtained not only for 
a zero longitudinal electric current (Michel 1973a; Mestel & 
Wang 1979), but also for a monopole magnetic field of a 
neutron star with a nonzero electric current (Michel 1973b). 
Another solution for a nonzero electric current was also dis- 
cussed (Beskin, Gurevich & Istomin 1983; Fitzpatrick & 
Mestel 1988ab; Lyubarskii 1990; Sulkanen & Lovelace 1990). 
As to the case of inclined rotator, the solution was obtained 
for the magnetosphere without a longitudinal electric cur- 
rent (Henriksen & Norton 1975; Beskin et al 1983; see also 
Mestel & Wang 1982). Further on, a more general approach 
including particle mass was discussed as well (Ardavan 1979; 
Bogovalov 1990, 1991; Mestel & Shibata 1994). 

However, the full analysis of the stream equation is far 
from being completed even for the force-free axisymmetric 
magnetosphere. The point is that the stream equation de- 
scribing the pulsar magnetosphere in the presence of a lon- 
gitudinal electric current becomes nonlinear. As a result, 
the comprehensive analysis of this equation meets certain 
problems (Michel 1991; Beskin, Gurevich & Istomin 1993). 
Nevertheless, in some cases the stream equation can be an- 
alyzed. 

In this work we consider the simplest approach in which 
the returning current fiows in a narrow layer along the sepa- 
ratrix, which is the boundary between open and closed field 
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line regions, and the current density in the open-line re- 
gion is constant. In this case, the stream equation describ- 
ing the magnetic field structure remains linear not only in a 
closed, but in an open magnetosphere as well. As a result, 
it is possible to write the general solutions describing the 
magnetic field structure in both regions, the main problem 
lying in matching these solutions along the separatrix. We 
shall show that the stability conditions on the separatrix 
allow us to obtain a unique solution of the problem. More- 
over, the hypothesis of the compatibility condition (Beskin 
et al 1983) connecting the potential drop in the pair cre- 
ation region and the longitudinal electric current flowing in 
the magnetosphere will be conflrmed. 



2 BASIC EQUATIONS 



Let us consider an axisymmetric stationary magnetosphere 
flUed with plasma. This means a zero longitudinal electric 
field in the magnetosphere 

(EB) = 0. 

In this case, the magnetic and electric fields can be repre- 
sented as 

B = ^l^ + ^e,, (1) 

VJ czu 

E^-^Vf, (2) 
c 

where c is the velocity of light, and Cl is the angular ve- 
locity of the star. Here f{'OJ, Z) is the fiux function, which 
is constant along magnetic field lines, and g{zu,Z) deter- 
mines the whole electric current I fiowing inside the tube 
/(ra7, Z) = const: 

21 

Finally, the value ^f(t^,Z) is the angular velocity of the 
plasma. Indeed, determining the drift velocity from equa- 
tions (111) and (PI) we have 



[E X B] 



where i\\ is a scalar function. Hence, the particle velocity is 
a sum of the corotation velocity and the slide velocity along 
magnetic field lines. Using now the Maxwell equation 

[V xE\= 0, 

which holds in the stationary case, one can see that the an- 
gular velocity Q.f must be constant on the magnetic surfaces 
/ = const 



Qf = !^f(/). 



(3) 



This is a well-known Ferraro izorotation low (Alfven & 
Folthammar 1963). It is clear that Qf can differ from the 
angular velocity O of the neutron star only on field lines 
passing through the polar cap region with a longitudinal 
electric field, i.e. through the particle creation region. As a 
result, in the open magnetosphere the angular velocity Qf 
must be smaller than that of the star. On the other hand, 
we have Qf = n in the region of closed field lines. 



Inserting now the electric and magnetic fields (|l]) 
and (^ into the force-free condition 

E -VE + [\/ X B] X B = 0, (4) 
we obtain from the ip component of equation (^) 
V/ X Vg = 0. 

Thus, the electric current g{vD,Z) must be constant on the 
magnetic surfaces as well 



9if)- 



(5) 



As to the other components of equation (j^, they just give 
us the stream equation describing the stable configuration 
of the magnetic field lines. In the dimensionless form 

X = Qzu/c; z = Q,Z/c; to = Q,f/Q 

we have 



A/(l-c^V) 



2 df Ag 2 dtj 2 „ la\ 

^a^-^dj + ^'d/^^^) =°- 

This is the general form of the force-free equation de- 
scribing a stable magnetic field configuration in the axisym- 
metric case (Okamoto 1974). We see that in equation (^ 
the two last terms depending on the longitudinal electric 
current <?(/) and angular velocity Q.F{f) can be nonlinear. 
On the other hand, in the closed region (region 1), where 
longitudinal electric currents are absent [g = 0), and the 
angular velocity is Q.f = f2, the equation becomes linear 
(Mestel 1973; Michel 1973a) 



A/W(l-x2)+2^/ 



(1) 



„ 0. (7) 

X ox 

Moreover, assuming that in the open field region (region 2) 
Qf = Q{1~Po), (8) 
9if) = iof, (9) 
where io and /3o are constants, we obtain the linear equation 



A/(2) [l_^2(^_^^^)2]_^ 2 9/(2) 



da 



(10) 



for the open region as well. Physically, the parameter /3o 
corresponds to the potential drop $ in the double layer near 
the star surface, i.e. in the particle creation region 



/30 = */*ma:c, 

where 



(11) 



(12) 



[R is the neutron star radius, Bq is the magnetic field on 
the star surface) is the maximum potential drop that can be 
realized near the star surface inside the open field line region 
(Ruderman & Sutherland 1975). For * ^ '^max, A ~ 1 
we have f^F ~ 0, i.e. there is no plasma rotation on the 
open field lines. Actually, the value of Po can be determined 
by the pair creation process in the double layer near the 
magnetic polar caps (Ruderman & Sutherland 1975; Arons 
1983; Gurevich & Istomin 1985). On the other hand, we have 



3GJ 



(13) 
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where 



3gj 



2-K 



is the Goldreich-Juhan current density. The condition ^ 
just means that the current density in the open field region 
is constant. Thus, with this approach we have two parame- 
ters io and /3o, the latter being actually determined from a 
certain pair creation mechanism near the star surface. The 
question arises whether it is possible to determine the elec- 
tric current io as well. 

We note, that the stream equation (^) is valid only in- 
side the light surface when \E\ < \B\. Otherwise, one can 
not neglect particle inertia in the force-free equation 
That is why, for example, the region of the closed-line mag- 
netosphere can not spread beyond the light cylinder x = 1. 
The closed-line region is limited by a zero point of the mag- 
netic field Xt, which is the point of intersection of the sepa- 
ratrix with the equatorial plane. Therefore, we must impose 
the condition a;, < 1 to construct closed-line magnetosphere. 
As for open-line magnetosphere, the light surface in this re- 
gion lies beyond the cylinder x — 1/(1 — /3o) because of the 
slowing-down of the magnetosphere rotation. 

Both equations (Q) and (|lo|) contain no coordinate z. As 
a result, the general solutions can be represented as (Mestel 
& Wang 1979; Beskin et al 1983; Mestel & Pryce 1992) 



/(3 



ip{\)R{x, X) cos XzdX. 



(14) 



As to the radial function R{x,X), it satisfy an ordinary 
second-order differential equation (see below). 

In the open region it is convenient to introduce new 
variables 

xi=x{l~l3o); zi = z{l-f3o). 

One can check that in these variables equation ( p^ only 
depends on one parameter 



(1 " /3o)2 



Moreover, as was shown by Beskin et al (1983, 1993), for an 
open field line region we must omit all harmonics between 
A = and X — a\^^ in the integration (p^). Otherwise, the 
magnetic field along the rotation axis will have sinusoidal 
oscillations, i.e. magnetic islands. As a result, the general 
solution in an open field region can be rewritten in the form 



f'Hx,z)^ 



(f2(X)R2{xi, X) ( 



(A^ -fell) 2 21 



dX, (15) 



where now the radial function _R2(2;i,A) is the solution of 
the ordinary differential equation 

d^7?2(cci,A) l + xf dR2{xi,X) 



dxl 



Xl(l 



ai- 



dxi 
-R2(2:i, A) 



(16) 



3 BOUNDARY CONDITIONS 



Let us now consider the boundary conditions for equa- 
tions (0) and (|l^). First of all, at small distances the so- 



lutions of both equations must be matched to the magnetic 
field of a neutron star. More exactly, the normal component 
of the magnetic field Bn must be continuous on the star 
surface (Bogovalov 1991). Hence 



f{x,z) = fstar{R,0); x = RcosO; z — RsinO. 
For example, for a dipole magnetic field we have 

f{x,z) 



(a;2 -f 22)3/2 ' 



(17) 



(18) 



where m is the magnetic dipole of a neutron star, so 

fstariR,e) = mshi^e/ R. 

It is also necessary to impose the regularity conditions 
on the singular surfaces of equations (^^ and (Q) 



2df^ 
X dx 



l3o 



df 



(1) 



da 



0; 



(19) 



(20) 



On the other hand, if the zero point is located inside the light 
cylinder, it is not necessary to impose the condition (^o|) for 
the closed-line region. 

Next we will impose an obvious symmetry condition for 
the magnetic field in the closed field line region 



df 



(1) 



dz 



0; 



0, 



(21) 



and the condition of zero magnetic flux on the z-axis for the 
open field line region 



/ 



(2) 



0; x^O, 



(22) 



The last two conditions are automatically fulfilled if we rep- 
resent the solution by (^) and v'2(A) does not have any 
singularities in the upper half-plane of the complex A. 

Finally, we must match the solutions in the open and 
closed regions. First of all, it means that the boundary of the 
closed region f^^\x,z) — /, must coincide with the position 
of the open region field line /'^^ (x, z) = f,, that is. 



r\x,z,{x)] = ,f'^[x,z,{. 



(23) 



where z = z^ix) is the position of the separatrix / = /, . 
One can see that our approach at this point differs from 
that of Lovelace & Sulkanen (1990) who proposed the gap 
between open and closed regions. Moreover, as pointed out 
by Okamoto (1974) and Lyubarskii (1990), for stability of 
the separatrix separating open and close regions it is neces- 
sary that the value of L = B^ — E"^ be continuous through 
this separatrix 

L^^\x,z,{x)] = L^^\x,z,{x)]. 

Using definitions (1), (2) and relations (8), (11), and (13), 
we obtain 

x^L = \yf^'^\''[l-x^] = 
\Vf^^^\'[l-x\l-Pof]+ilfl 



(24) 



This condition results from the integration of the full stream 
equation (^) over the narrow boundary f ~ ft- This is 
the consequence of nonlinearity of the general stream equa- 
tion (^. As a result, to construct the solution it is necessary 
to specify two parameters io, /3o, and the condition ( p^ on 
the star surface. All the other characteristics, in particular. 
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the position and the structure of the zero point can be de- 
termined from the above boundary conditions. 

At the very end, it is necessary to stress that formally, to 
fulfil the condition (2^) at the zero point Xt, we can expect 
X, < 1, (a/(^V9a;)|^^ / 0, and {d'' f^""^ /d^z)\^^^ = 0. 

In other words, here | V/*-"*^' |^ [1 — x'i] — i^f^ ( Lyubarskii 
1990). Therefore the angle 2i9 between intersecting separa- 
trices will equal 2f? = n. However, there is always a narrow 
transition layer 5f along the separatrix. As a result, the 
magnetic structure in this region remains X-point as for a 
magnetosphere with io = and /3o = 0. The only difference 
is that, as one can easily see from (^), the value A/|^^^ 
is equal to zero if x, < 1, so 2i9 = 7r/2. Indeed, supposing 
that the returning current density is constant in this layer, 
we can write <?(/), ^{f) and L for /* — Sf < f < f* as 

g = ioif,~5f){f,-f)/5f, 
L^ = l-/3o(/*-/)/<5/, 




0.00 0.04 



Po 



Figure 1. Forbidden and permitted regions in the io — A) dia- 
gram. The solid curve corresponds to the compatibility condition 
*o(/3o)i while the dashed curve represents the dependence (26). 



2 r 

x L 



iv/i= 



-A) 



5f 



+ 



ilU,-5fY 



Hence, L is continuous through the transition layer and the 
separatrix. On the other hand, if 5f tends to zero, as is con- 
sidered in this work, then we have to take into account that 
the condition ( p^ is not valid near to zero point. But this 
is natural, because the stream equation becomes essentially 
two-dimensional near the zero point. 



FORMULATION OF THE PROBLEM AND 
NUMERICAL RESULTS 



In their previous paper Beskin et al (1983) assumed that de- 
spite the presence of an electric current io / in the region 
of open field lines the structure of closed magnetosphere re- 
mains the same as for io = 0, /3o = 0. In particular, the zero 
point is proposed to lie on the light cylinder a; = 1. It was 
shown that for these assumptions the values jq and /3o must 
be bounded by the compatibility condition 



Po = /?,(jo), 
where 

P,{io) = 1 - 



1/2 



(25) 



(26) 



and imax ~ 1.58. It means that separatrix field lines f = f* 
in open and closed regions coinside if the relation (26) is 
carried out only. The dependence (^) is shown in Fig. 1 by 
the dashed line. The relation ( |26| ) can be derived directly 
from the stream equation (p^). Indeed, assuming that the 
field line / = /* passes near the zero point a;* = 1 (where 
df^'^^ /dx = 0), we immediately obtain the condition ( |2^ ) 
with t^ax = (A/|j;=j. )//*■ The goal of the present paper 
is actually to check the compatibility condition ( |25[ ) in a 
more self-consistent approach, i.e. to include the disturbance 
of the closed region and the stability condition (p^). 



As has already been shown, equations (^, ( p^ ) and the 
boundary conditions (|l7[), (p^-(|2^ give a solution of the 
problem if it exists. On the other hand, we do not know 
beforehand the position of the separatrix. It can be deter- 
mined from the solution. In other words, the function <^2(A) 
in (p^, which determines, in particular, the position of the 
separatrix and the zero point, is unknown. However, the 
magnetic field has a dipole asymptotic behavior when x, z 
tend to zero. For the representation ( p^ ) and the star mag- 
netic dipole m = 1, this assumption corresponds to 



lim 7i2 (2:1 , A) — —\xKi{\x), 

x->-0 TV 

because we have for a dipole field (^ 

x^ 2 



(27) 



(2:2+^2)3/2 

Hence, 



lim V'2(A) 

A — *oo 



XxKi {Xx) cos XzdX. 



(28) 



On the other hand, from (|l9|), ( |l5|) we have for the function 
R2{xi,X) 



2 dR2 

xi dxi 



oiiR-z 



0; 



Xl 



(29) 



As a result, for given parameters io and /3o the func- 
tion 1^2 (A) completely determines the open-line magneto- 
sphere. Then using equation (Q) and the boundary condi- 
tion ([23[), ( |2l| ) we can determine the closed-line magneto- 
sphere as well. 

In the present work, for any pair of parameters from 
the io-Po plane we are looking for the best function 952 (A) 
to satisfy the condition (P4). In other words, to obtain the 
continuity of L = — E through the separatrix, for any 
io, /3o we minimize the functional 



S[^p2{X)]= max \l'^'^\x, z,{x)] 



L'^'\x,z,{x)]\, (30) 



where 5x takes into account that the condition ( p4[ ) is not 
valid in the vicinity of the zero point. In more detail, we 
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expand the function tp2{\) in a finite functional series 



tf2 



(A) = 1 + ^C„(^„(A) 



to reduce tlie problem of functional minimization to the 
problem of minimization of function with many variables, 
which are coefficients of the expansion. In our work we use 
iV = 2 and 



<^i(A) = 



<f>2{X) 



1, AG [0.2, 0.6] 

0, A ^[0.2, 0.6] 

1, AG [1.5, 3.0] 
0, A ^[1.5, 3.0] 



So, the function </52(A) is different from one on two inter- 
vals A G [0.2,0.6] and A G [1.5,3.0]. The resuh does not 
sufficiently depend on what intervals to choose. The func- 
tion <j)\{\), that is nonzero for smaller A-interval, and the 
function (^2 (A), that is nonzero for greater A-interval, allow 
us to minimize the functional (|o|) for x ^ 1 and j: ^ 1 
respectively. 

As has already been said, the zero point must lie inside 
the light cylinder x, < 1. On the contrary, the open region 
can not be matched to the closed one because the closed- 
line magnetic field may not spread beyond the light cylinder 
a; = 1. As a result, a self-consistent solution of the problem 
can be obtained not for all the parameters io f3o- This is 
illustrated in the io-Po diagram in Fig. 1. The solid line is 
the boundary of the forbidden region of the parameters io, 
f3o. For the forbidden region, the zero point of the magnetic 
field moves outside the light cylinder. 

As for the permitted region, numerical calculations 
show that for a given /3o we get the best conformity to the 
condition (^) if we take io that lies just on the boundary 
between the forbidden and the permitted regions. Simulta- 
neously, the zero point of the magnetic field just lies on the 
light cylinder x, = 1. Moreover, the minimum of the elec- 
tromagnetic field energy 



x'^{l + x'^)\Vf'-^''fdzdx + 



|V/(^)|^(l + x^(l-/3o)^)+zg(/(^))^ 



dzdx 



z,{x,) 

corresponds to the boundary curve as well. To illustrate the 
above statements we show the dependence of the zero point 
Xt, the value of the functional S[if2{X)\ (30) corresponded 
to the drop of L = — E"^ and the electromagnetic field 
energy Eem from io for the /3o = 0.05 in Figs. 2a, 2b and 2c 
respectively. As is seen from Fig. 2b, we can not find the 
exact function (^2 (A) to obtain zero drop of L because we 
would have to expand (^2 (A) in an infinite functional series 
in this case. 

Thus, we can conclude that the compatibility condition 
/9o = /3*(io) does actually exist. It is represented by the 
boundary solid line in Fig. 1. One can see that the difference 
between the dashed and solid lines is very slight for small 
values of io, /3o, and hence we conform the dependence ( |2^ ) 
for a small current density and a potential drop. 



X. 

1,00 
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0,97 
0,96 
095 



0,5 0,6 0.7 0.8 0.9 1.[ 



s[(pa): 



_l I I 1 I I L. 



0.7 0.8 0.9 1.0 



6700 - 



Figure 2. The dependence of the zero point x, (A), the value 
of the functional 5(1^2 (A)] (30) corresponded to the drop of L = 
— (B) and the electromagnetic field energy E^m (C) from 
io for the /3o = 0.05. One can see that for io = 0.52, which 
corresponds to the compatibility condition, we have x, = 1 and 
the minimum of the values S[(p2{\)], Eem- 



The structure of the magnetosphere for /3o — 0.05 and 
io = 0.52, which corresponds to the compatibility condi- 
tion, is shown in Fig. 3. We see that the zero point of the 
magnetic field lies just on the light cylinder x = 1 and the 
angle 2i? between the intersecting separatrices corresponds 
to tan 1? = 1/2. Therefore, we confirm the hypothesis by 
Beskin et al (1983) that the shape of the closed region in 
the presence of all electric current and a potential drop does 
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Figure 3. The structure of the axisymmetric magnetosphere of a 
neutron star for the potential drop and the current density which 
are connected by the compatibility condition. The zero point is 
located on the light cylinder x = I and the angle 2i? between the 
intersecting separatrices corresponds to tan^ i? = 1/2. 





|3„=0.05 




i,=0.78 








1 1 


1 1 1 



Figure 4. The structure of the axisymmetric magnetosphere of a 
neutron star for the potential drop and the current density which 
are far away from the compatibility condition. The zero point is 
located inside the light cylinder x = 1. 



not differ greatly from that for the zero current. On the 
other hand, in Fig. 4 we show an example, where the zero 
point lies inside the light cylinder {x, < 1). This can be real- 
ized for the parameters io and /3o far from the compatibility 
condition. 



5 DISCUSSION 



Thus, we have shown that the compatibility condition be- 
tween current and potential drop remains true even if we 
include the boundary condition (^) into consideration. On 
the other hand, it was shown by Beskin et al (1983) that the 
relation (^) remains, in general, the same for an inclined 
rotator as well. For example, the compatibility relation has 
now the form 



■%ax (x) 



1/2 



(31) 



where f3max and imax depend now on the inclination angle X- 
As a result, the longitudinal current flowing in the neu- 
tron star magnetosphere is not a free parameter and deter- 
mined by the particle creation mechanism in the polar re- 
gions. On the other hand, it is the longitudinal current that 
causes deceleration of pulsar rotation. For example, the total 
energy losses in the model proposed by Beskin et al (1983) 
is 



Irftn = Wtot 



-ioiPo) cosx, 



(32) 



where x is the angle of pulsar inclination, Ir is a moment of 
inertia of the pulsar and 1.59 < /*(x) < 1-95. 

The full analysis of observational data goes beyond the 
scope of this paper. Ifowever, we briefly recall the main 
predictions which result from the existence of the Ohm's 
Law (H). 

1. As we understand the nature of pulsars, their radio 
emission results from a secondary plasma which is produced 
by a longitudinal electric field near the pulsar surface. There- 
fore, the condition 4'(P, B) = "ifmax just determines the 
maximum period of radio pulsars. For example, in the model 
with large enough work function (Ruderman & Sutherlend 
1975; Gurevich & Istomin 1985) we have Pmax ~ Is, which 
is in agreement with observations. Using the compatibility 
condition io = io{/3o) and (|ll|), (p^, (p2[), we can rewrite 
the inequality P < Pmax in the form Q < 1, where 

.p. 11/10/ p X-4/10 



The parameter Q ~ io, which can be determined from ob- 
servations, is very convenient for expressing the basic pul- 
sar characteristics (Beskin et al 1984; Taylor & Stinebring 
1986; Rankin 1990). For example, in the hollow cone model 
the ratio of the inner radius Tin and the height H of plasma 
generation region to the polar cap radius Ro of a pulsar is 

rin ^7/9. H 

-Ro^"^ ' -Ro^^- 

Thus, we can conclude that pulsars with Q ~ 1 have a thin 
cone of emission and hence they have a two-peak radio emis- 
sion profile. It is these pulsars that have emission nonregu- 
larity such as nulling, mode switching, etc. On the contrary, 
pulsars with Q have a stable one-peak profile emission. 
Such situation is just corresponds to observations (Taylor & 
Stinebring 1986; Beskin et al 1993). 

2. One of the main radio pulsar parameters that char- 
acterizes pulsar rotation deceleration is a braking index 
ni,r = f2fi/ which is directly available from observations 
(Lyne & Graham-Smith 1990). Unfortunately, the braking 
index is only known for a few radio pulsars. In particu- 
lar, = 2.24 for PSR 0540 - 693 and = 2.84 for 
PSR 1509 — 58. In the case of vacuum magnetodipole en- 
ergy losses we have Ubr =3-1- 2tan~^Xi which can not 
explain the observations. On the other hand, for the cur- 
rent losses ( p^ ) and the compatibility relation jo(/3o) ( |26| ) 
we have n^,r = 1.93 -1-1.5 tan'^ Xi which is in good agreement 
with observational data. 

3. In the Ruderman-Sutherlend model (1975) the longi- 
tudinal current is equal to the Goldreich- Julian one. Hence, 
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the energy of particles hitting against the pulsar polar cap is 
high enough to cause an intensive X-ray emission. This con- 
tradicts observational data. Nevertheless, the X-ray emis- 
sion is only observed to fast-period pulsars which have a 
considerably smaller current. Then, in accordance with the 
model discussed and the compatibility condition, the energy 
of pulsar cap heating is Wx ~ Q^Wtot (Beskin et al 1993), 
which is significantly less then the generally accepted esti- 
mates because for these pulsars sa = (j/jaj)^ <S 1. It 
is necessary to stress that in the Arons model (1993) with a 
small work function this difficulty is absent. 

4. If a pulsar energy lost is caused by longitudinal cur- 
rent, then the braking angular momentum K is opposite to 
the magnetic dipole of a neutron star. As a result, the value 
sin X remains constant during radio-pulsar evolution. Con- 
sequently the angle x between dipole axis and rotation axis 
of the pulsar is increasing during pulsar life, while for mag- 
netodipole energy lost it is decreasing. 

Thus, as we see, the main predictions of the theory 
with a large work function and the compatibility condi- 
tion /9o = /9*(*o) at least do not contradict observational 
data. Moreover, to confirm the theory, there have recently 
appeared some indirect observational results such as the ab- 
sence of evolution of radio-pulsar magnetic field and the sta- 
tistical conclusion that the initial pulsar periods are about 
0.1 — 1 s rather than 1 — 10 ms as for Crab and Vela young 
pulsars. Of course, the direct measurement of the evolution 
of pulsar inclination angle x will be the key-experiment. Un- 
fortunately, such an experiment is impossible to carry out 
at the present time. 
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